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Abstract. We use the open clusters (OCs) with known parameters available in the WEBDA database and in 
recently published papers to derive properties related to the disk structure such as the thin-disk scale height, 

■ displacement of the Sun above the Galactic plane, scale length and the OC age-distribution function. The sample 
totals 654 OCs, consisting basically of Trumpler types I to III clusters whose spatial distribution traces out the 

j. , local geometry of the Galaxy. We find that the population of OCs with ages younger than 200 Myr distributes in 

■ the disk following an exponential-decay profile with a scale height of Zh = 48±3 pc. For the clusters with ages in the 
range 200 Myr to 1 Gyr we derive Zh = 150 ± 27 pc. Clusters older than 1 Gyr distribute nearly uniformly in height 

„ , from the plane so that no scale height can be derived from exponential fits. Considering clusters of all ages we 

^ ' obtain an average scale height of Zh — 57±3pc. We confirm previous results that Zh increases with Galactocentric 

distance. The scale height implied by the OCs younger than 1 Gyr outside the Solar circle is a factor ~ 1.4 — 2 
larger than Zh of those interior to the Solar circle. We derive the displacement of the Sun above the Galactic plane 
as zq = 14.8±2.4 pc, which agrees with previous determinations using stars. As a consequence of the completeness 
effects, the observed radial distribution of OCs with respect to Galactocentric distance does not follow the expected 
exponential profile, instead it falls off both for regions external to the Solar circle and more sharply towards the 
Galactic center. We simulate the effects of completeness assuming that the observed distribution of the number 
of OCs with a given number of stars above the background, measured in a restricted zone outside the Solar 
«4 ' circle, is representative of the intrinsic distribution of OCs throughout the Galaxy. Two simulation models are 

considered in which the intrinsic number of observable stars are distributed: (i) assuming the actual positions of 
the OCs in the sample, and (ii) random selection of OC positions. As a result we derive completeness-corrected 
' radial distributions which agree with exponential disks throughout the observed Galactocentric distance range 

\ 5-14 kpc, with scale lengths in the range Rd = 1.5 — 1.9 kpc, smaller than those inferred by means of stars. In 

particular we retrieve the expected exponential-disk radial profile for the highly depleted regions internal to the 
Solar circle. The smaller values of Rd may reflect intrinsic differences in the spatial distributions of OCs and 
stars. We derive a number-density of Solar-neighbourhood (with distances from the Sun d© < 1.3 kpc) OCs of 
Pq = 795 ± 70 kpc -3 , which implies a total number of (Trumpler types I to III) OCs of ~ 730 of which ~ 47% 
would already have been observed. Extrapolation of the completeness-corrected radial distributions down to the 
Galactic center indicates a total number of OCs in the range (1.8 - 3.7) x 10 5 . These estimates are upper-limits 
because they do not take into account depletion in the number of OCs by dynamical effects in the inner parts of 
the Galaxy. The observed and completeness-corrected age-distributions of the OCs can be fitted by a combination 
of two exponential-decay profiles which can be identified with the young and old OC populations, characterized 
by age scales of ~ 100 Myr and ~ 1.9 Gyr, respectively. This rules out evolutionary scenarios based on constant 
star-formation and OC-disruption rates. Comparing the number of observed embedded clusters and candidates in 
the literature with the expected fraction of very young OCs, derived from the observed age-distribution function, 
we estimate that 3.4-8% of the embedded clusters do actually emerge from the parent molecular clouds as OCs. 

Key words. (Galaxy:) open clusters and associations: general; Galaxy: structure 
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1. Introduction 

Understanding the true shape of the Milky Way is a chal- 
lenging task which depends essentially on the availability 
of as complete as possible Galactic surveys of the spa- 
Send offprint requests to: C. Bonatto 
Correspondence to: charles@if.ufrgs.br 



tial distribution of stars, star clusters, Hi gas, molecular 
clouds, etc. These observations can be used to better con- 
strain theoretical mass models of the Galaxy. The modern 
view of the Milky Way pictures it as composed basically of 
the center and nuclear bulge, a thin disk hosting most of 
the OCs, a thick disk containing the old OCs, and the ex- 
tended halo hosting the old stars, white dwarfs and glob- 



2 



C. Bonatto et al.: Spatial distribution of open clusters 



ular clusters (e.g. Majewski [1993 ). The gas and stellar 
density of the disk is usually expressed as a combination 
of exponential-decay profiles for the horizontal and ver- 
tical directions, p(r,z) oc e~( r / RD ' e~^ z ^ Zh \ where Rd is 
the scale length and z h the scale height (e.g. Binney & 
Tremaine HWj) . 

The kinematics and spatial distribution of large num- 
bers of stars have been the best used probes of the 
Galaxy structure. For instance, M dwarfs are found to 
distribute vertically with relatively large scale heights, 
Zh ~ 300 pc, in contrast to the younger A-type stars 
with Zh ~ 100 pc (Faber et al. I1976j) . Giants and sub- 
giants have scale heights typically between 200 and 500 pc, 
though the uncertainty is larger than that of the dwarf lu- 
minosity class (e.g. Bahcall & Soneira 1981). Using star 
counts from the Sloan Digital Sky Survey taken from 
two high-latitude (49° < |b| < 64°) samples north and 
south of the Galactic plane Chen, Stoughton, Smith et 
al. lj2"nt]l"|l derived z h = 330 ± 3 pc for the old thin disk 
and Zh = 580 — 750 pc for the thick disk. They suggest 
that the thick disk formed through heating of a preexist- 
ing thin disk caused by the merging of a satellite galaxy. 
Based on the kinematics of sdB stars de Boer et al. I|1997fl 
and Altmann et al. (2004) found that these stars are char- 
acterized by two distributions, a disk one with Zh ~ 1 kpc 
and an extended one with Zh ~ 7 kpc. In recent years con- 
siderable amounts of data on parallax and proper motion 
of stars has allowed to probe in detail the Galactic struc- 
ture in distances of up to 2 kpc from the Sun. In particular 
Kaempf, de Boer & Altmann (2005} using the kinemat- 
ics and spatial distribution of red horizontal-branch stars 
selected from the HIPPARCOS catalogue derived a disk 
scale height of Zh ~ 0.6 kpc. 

With respect to the vertical position of the Sun (z ), 
asymmetries in star counts above and below the Galactic 
plane have put it displaced a few parsecs to the North. 
Several methods used in the past constrained the value of 
z© to the range 4-40 pc (see, e.g. the account in Cohen 
1995|. More recently, different sets of data have been used 
to obtain more accurate values of z©. Cohen (1995), us- 
ing IRAS point-source counts at 12 and 25 /im obtained 
z Q = 15.5 ± 0.7 pc. Hammersley et al. fl_9_9_5J examined 
star counts and surface brightness maps from COBE, 
IRAS and the Two-micron Galactic Survey to derive 
z = 15.5 ± 3 pc. The high-Galactic latitude data of Chen, 
Stoughton, Smith et al. ( 20Q]J resulted in z = 27±4pc. 

Measurements of the horizontal scale length by means 
of stars place Rd in the range from 2.2 to 3.5 kpc (de 
Vaucouleurs & Pence 119781 Knapp, Tremaine & Gunn 

nwm McCuBkev iissfl . 

Open clusters are formed in, and distribute through- 
out the disk. Interactions with the disk and the tidal pull 
of the Galactic center/bulge tend to destroy the poorly- 
populated OCs in a time-scale of a few 10 8 Myr (Bergond, 
Leon & Guibert 120 OT)) . The dynamical disruption is more 
critical particularly for the OCs more centrally located. 
The survivors, however, may reach large vertical distances 
from the plane spanning the thick disk. Although the num- 



ber of known OCs is small compared to stars, it is rela- 
tively simple and accurate to derive distance and age for 
these objects. In this sense, OCs - in particular their spa- 
tial distribution - may serve as a direct probe of the disk 
structure. By means of a sample of 72 OCs Janes & Phelps 
<|r$9Hl found that the young OCs (age < 800 Myr) are dis- 
tributed almost symmetrically about the Sun with a scale 
height perpendicular to the Galactic plane of Zh « 55 pc, 
while for the old OCs they found Zh « 375 pc. 

It is interesting to point out that since Janes & Adlcr 
(1982) study the number of OCs with derived parame- 
ters increased nearly 50%, taking as reference the clusters 
currently with parameters in the WEBDA 1 OC database 
(Mermilliod 1996). However, in the recent revision of OCs 
Dias et al. (2002) report a total of 1537 catalogued OCs. 
This clearly indicates the need to explore more systemat- 
ically the available OC databases. WEBDA gathers data 
from a variety of sources, which implies different observa- 
tional instrumentation and analysis methods and, conse- 
quently uncertainties in the values of parameters. We deal 
with the non- uniformity of the parameters by (a) work- 
ing with bins wide enough to take into account most of 
the uncertainties, and (b) assuming la Poisson errors to 
assign statistical significances to the results. 

In this paper we build an OC sample composed of the 
WEBDA objects added to 10 other OCs studied in recent 
papers. We intend to use this OC sample, particularly the 
z-measurements, to derive the disk scale height (zh), the 
displacement of the Sun above the Galactic plane (zq), 
and the horizontal scale length. In addition, we use the 
present OC sample to study statistical properties related 
to the reddening, distance from the Sun, location in the 
Galaxy, completeness and, particularly, to derive the in- 
trinsic OC age-distribution function. 

Previous statistical works on the relation of OC pa- 
rameters with Galactic structure are, e.g. Lynga (1982), 
Janes & Phelps fiMfy . Nilakshi et al. QMty . and Tadross 
et al. I|2002[) . The main results of Lynga (1982), related 
to the present work, are (i) a correlation of cluster age 
with Galactocentric distance, (ii) an exponential OC dis- 
integration time scale of ~ 10 s yr, (Hi) displacement of 
the Sun over the Galactic plane of 20 pc, and (iv) a 
z-distribution which differs only for the oldest clusters. 
Nilakshi et al. I|2002)) found that OCs with Galactocentric 
distances larger than 9.5 kpc have larger sizes. Tadross et 
al. ( 2002) found (i) evidence that older clusters are found 
preferentially in the outer parts of the disk, and (ii) that 
most of the OCs younger than 32 Myr are concentrated on 
the Perseus arm. 

This paper is organized as follows. In Sect.Elwe present 
the OC sample and discuss statistical properties related 
to the spatial and age distributions. In Sect. |3] we dis- 
cuss qualitative aspects affecting sample completeness. In 
Sect. 0]we derive the displacement of the Sun above the 
Galactic plane and the observed scale height. In Sect. El we 
simulate the effects of sample completeness. In Sect.Elwe 

1 http://obswww.unige.ch/webda 
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derive completeness-corrected radial and vertical distri- 
bution functions, and discuss properties of completeness- 
corrected age distribution functions. In Sect. [7| inferences 
on the total number of OCs in the Galaxy are presented. 
Concluding remarks are in Sect. [SJ 

2. Statistical characterization of the OC sample 

The OCs used in the present study are those with Galactic 
coordinates, distance from the Sun, reddening and age cat- 
alogued in WEBDA. To these we add 10 recently studied 
OCs (2 in Bica, Bonatto & Dutra l2T)0"3l 3 in Bonatto, Bica 
& Dutra UHHl and 5 in Ortolani, Bica & Barbuv r27)TI5| . 
The sample totals 654 OCs with known distance from the 
Sun (of which 636 have age determination). 

The significant number of OCs in this sample presents 
an opportunity to study properties of clusters in differ- 
ent age ranges separately. We selected the age ranges (i) 
age < 200 Myr (young OC population), (ii) 200 Myr < 
age < 1 Gyr (moderate-age OCs), and (Hi) age > 1 Gyr 
(old OCs) both for statistical and population representa- 
tivity purposes. The number of clusters in each age range 
is 402 (young), 148 (moderate-age) and 86 (old). 

2.1. Spatial distribution 

In Fig. ^ we show the spatial distribution of the OCs in 
the present sample projected onto the X and Y plane, 
separated according to age groups. Throughout this pa- 
per we use 8.0 kpc (Reid 1993) as the distance of the Sun 
to the Galactic center. In this figure conspicuous struc- 
tures of the local disk are apparent. For clarity we show 
the Solar circle in all panels. Note that in panels (a) and 
(b) the threshold between young and moderate ages has 
been changed to 60 Myr in order to better isolate OCs 
related to spiral arms. OCs younger than 60 Myr (panel 
(a)) trace out the local distribution of the Orion Arm 
(Georgelin & Georgelin[l970), in which the Sun is located 
close to the inner border. Another conspicuous feature de- 
fined by the young OCs is the gap towards the Galactic 
center at dec ~ 7.5 kpc related to the dip before the Sgr- 
Car Arm. Clusters older than 1 Gyr (panel (c)) are uni- 
formly distributed but with an asymmetry in the number 
of OCs towards the anti-center, a region where they are 
well-known to populate (Friel 1995). The all-ages sample 
is shown in panel (d) where traces of the Orion Arm, the 
gap at dec ~ 7.5 kpc and the radial asymmetry can still 
be detected. 

2.2. Observed age histogram 

The observed age histogram of the present OCs is shown 
in Fig. [21 with bins of 200 Myr which are wide enough to ac- 
commodate most of the uncertainties related to the differ- 
ent age-determination techniques included in the WEBDA 
sample. Two conspicuous peaks occur at ~ 100 Myr (char- 
acterizing the young OCs) and ~ 700 Myr (moderate-age 
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Fig. 1. Spatial distribution of the OCs, according to age 
groups. The Solar circle is shown by the dashed line, and 
the Galactic center is located at (0,0). Shaded area in 
panel (d): region encompassed by the restricted sample 
(Sect.[7J). For comparison purposes the coordinate defini- 
tions follow Janes & Phelps l|1994|l . 

OCs). A less-pronounced peak occurs at ~ 3.1 Gyr, which 
can be associated to the old OC distribution. 

The age-distribution function will be presented and 
discussed in Sect. 16.51 

2.3. Galactic coordinates 

In the left panels of Fig. |3|we examine the distribution of 
OCs in terms of Galactic longitude. In all age ranges the 
number of OCs drops off towards the Galactic center and 
increases roughly towards the third quadrant. These facts 
might reflect both an observational limitation in the sense 
that low-contrast clusters are more difficult to detect in 
directions projected against the Galaxy center/bulge (or 
more simply, observers tended to avoid crowded central 
fields - Sect.EJ) and the enhanced tidal-disruption proba- 
bility of OCs closer to the Galactic center/bulge (Bergond, 
Leon & Guibert HEDl . 

The distribution of OCs in terms of Galactic latitude 
is in the right panels of Fig. [21 As expected, most of the 
OCs are tightly concentrated close to the Galactic plane 
(panel (e)). However, the distributions become wider for 
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Fig. 2. Age histogram of the OC sample for age < 
3800 Myr. Peaks occur at ~ 100 Myr, ~ 700 Myr, and 
~ 3.1 Gyr. The inset shows the histogram for age > 
2800 Myr. 



older age groups (panels (f) to (h)). The distributions can 
be analytically described by the exponential-decay func- 
tion N oc e~l b / Zb l. We find z b = 2.10 ±0.20° for the young 
OCs (panel (f)), Zb = 3.47 ± 0.51° for the moderate-age 
OCs (panel (g)), and z b = 8.27 ± 1.13° for the old OCs 
(panel (h)). The all- ages distribution (panel (e)) is char- 
acterized by Zb = 3.14 ± 0.21°. The asymmetry observed 
in the old-age distribution with respect to Galactic lati- 
tude (panel (h)) may be accounted for by completeness 
(Sects. 01 and EJ- Because these clusters are not as concen- 
trated near the Galactic plane as the young OCs (Fig.0), 
the old OCs located above the Galactic plane become eas- 
ier to detect than the ones in the opposite side. A similar, 
but less-conspicuous effect can be seen in the distribution 
of the moderate-age OCs (panel (g)). Alternatively, this 
asymmetry might reflect an intrinsic difference in the dis- 
tribution of old OCs above and below the plane. 

2.4. Distance from the Sun and reddening 

In the left panels of Fig. 0] we show the distribution 
of OCs with respect to the distance from the Sun, in 
bins of 0.2 kpc. The young (panel (b)) and moderate- 
age (panel (c)) OCs present similar distributions with 
peaks at d© « 1.7 kpc and d Q « 1.3 kpc, respectively. 
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Fig. 3. Distribution of OCs in terms of Galactic longitude 
(left panels) and latitude (right panels). Solid lines: least- 
squares fit with the function N oc e~' b / Zb '. Panels (a) and 
(e): OCs of all ages. Panels (b) and (f): OCs younger than 
200 Myr. Panels (c) and (g): OCs with age in the range 
200 Myr - 1 Gyr. Panels (d) and (h): OCs older than 1 Gyr. 



Observational limitations probably affect the detection of 
low-contrast OCs for d Q > 1.7 kpc. The old OCs (panel 
(d)) are more evenly distributed (Fig.l, panel (d)). 

With respect to the reddening (right panels of Fig. 
the distribution of the young OCs (panel (f)) is wider than 
those of the moderate-age (panel (g)) and old (panel (h)) 
OCs. This fact probably can be accounted for by more 
significant internal reddening affecting young OCs. In ad- 
dition, excess interstellar reddening also affects more sig- 
nificantly young OCs, since these clusters are more tightly 
concentrated to the Galactic plane (panel (f) in Fig. 0J). 
The strong dependence of reddening on OC location in 
the Galaxy can be appreciated in panels (a) and (b) of 
Fig. El As expected, the reddening increases markedly to- 
wards the Galactic center/bulge (1st and 4th quadrants 
in panel (a)) and plane (the conspicuous peak in panel 
(b)). E(B — V) increases nearly linearly with the distance 
from the Sun up to rs 3 kpc (panel (c)); the scatter in- 
creases significantly for larger distances. The excess red- 
dening observed in the young OCs is apparent as well in 
panel (d). 
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Fig. 4. Distribution of OCs in terms of distance from the 
Sun (left panels) and reddening (right panels) . Panel ages 
as in Fig- EI 

2.5. Galactocentric distance 

In Fig. we examine histograms with the number of ob- 
served OCs in bins of Galactocentric distance. Similarly 
to the distance from the Sun, the old OCs (panel (d)) are 
more evenly distributed than the young (panel (b)) and 
moderate-age (panel (c)) OCs. The radial distribution of 
the all-ages OCs (panel (a)) is not symmetrical about the 
Galactocentric distance of the Sun, instead it presents a 
maximum at the Solar position and drops off more sharply 
towards the Galactic center than outside the Solar circle. 
OCs younger than 200 Myr and those with age in the range 
200 Myr - 1 Gyr present similar radial distributions. The 
observed radial distribution of the OCs clearly does not 
resemble that expected of an exponential disk (e.g. Binney 
& Tremaine ll987fl . We will return to this point in Sect.0 

3. Qualitative considerations on sample 
completeness 

In general terms, cluster detection-rates depend on fac- 
tors such as the stellar content (both in absolute numbers 
and the presence of giants), apparent magnitude at the 
turnoff, distance from the Sun, stellar background pro- 
jected onto the OC direction, etc. External factors such as 
systematic surveys with deep photometry are necessary as 
well to detect as many OCs as possible. In addition, the 
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Fig. 5. Dependence of reddening on Galactic longitude 
(panel (a)), latitude (panel (b)), distance from the Sun 
(panel (c)) and cluster age (panel (d)). Each point repre- 
sents an OC. 



identification of a stellar concentration as a star cluster 
depends on the analysis of the colour-magnitude diagram 
morphology, the structure of radial-density profiles, and 
proper-motion data, when available (e.g. Bonatto & Bica 
2005, Bica & Bonatto EUUSal Bica & Bonatto I2005EI and 
references therein). 

The present sample is not complete, particularly with 
respect to OCs at large distances from the Sun and/or low- 
contrast clusters. The cluster data in WEBDA come from 
photographic, photoelectric or CCD photometry taken 
with different instrumentation over the last 50 years. This 
heterogeneity of data sources precludes objective infer- 
ences on the faint magnitude-limit of the stars, total num- 
ber of stars, surface brightness or density contrast in the 
present OCs. Basically the WEBDA sample contains OCs 
of the denser Trumpler types I to III (Ruprecht 1966), be- 
cause observers naturally directed their attention to this 
kind of object. 

The following argument illustrates the difficulties as- 
sociated with the identification of low-contrast clusters in 
different regions of the Galaxy. Assume an OC whose ra- 
dial density profile can be described by the two-parameter 
King l|1966[) law, er(r) = cr bg + \a ; where Cbg is the 

background surface density of stars, ctok is the central den- 
sity of stars (above the background), and R c is the core 
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Fig. 6. Distribution of OCs in terms of Galactocentric dis- 
tance. Panel ages as in Fig. OH 



radius. In fact, the background represents the contamina- 
tion of non-cluster stars projected onto the OC angular 
area, integrated along the direction towards the OC, fore 
and background stars alike. In this sense the background 
density of an OC is a function of I and b only, and not of 
distance from the Sun (Sect.EJ. 

We define the contrast parameter as S c = — 
1 + Integration of <r(r) from the center to the lim- 

iting radius (Rn m ) and considering that Rii m > 5 x R c 
(Bonatto & Bica 2 005[) . the total number of cluster stars 
(background-subtracted) is N* w 27rer KRc l n (Riim/R. c )- 
Thus, the contrast parameter as a function of the back- 
ground density and effective number of stars can be ex- 
pressed as 



1 



(N * /a bg 



27rln(R lim /R c )R2 



(1) 



Extractions of 2MASS 2 photometry in directions 
which do not intercept detected clusters show that the av- 
erage background stellar density in low-latitude (\b\ < 5°) 
regions towards the Galactic center/bulge is a factor ~ 10 
as large as that towards the anti-center (Sect. EJ. Thus, 

2 The Two Micron All Sky Survey, All 
Sky data release (Skrutskie et al. ^997), 

http://wwwApac.caltech.edu/2mass/releases/allsky/ 



an anti-center cluster with 8 C ~ 11 would have S c ~ 2 
when projected against the center/bulge. In this sense, if 
the typical detection limit of WEBDA OCs is restricted 
to, say S c > 2, surveys towards the Galactic center/bulge 
would miss the whole class of OCs with 8 C in the range 
2 to 11 which in principle would be detected towards the 
anti-center. 

To give an idea of the values S c takes on with actual 
OCs wc quote here some examples taken from Bonatto 
& Bica (2005) and Bica & Bonatto l|2005b|l . derived with 
2MASS data. The high-Galactic latitude OCs NGC 188 
(b = +22.39°; d ~ 1.7 kpc; M - 3800 M Q ) and M67 
(6 = +31.89°; d - 0.9 kpc; M - 990 M Q ) present 5 C = 
22.5 ± 2.9 and 5 C — 34.7 ± 5.7, respectively, while the 
populous and rather low-latitude cluster NGC 2477 (b = 
-5.82°; d© ~ 1.2 kpc; M ~ 5300 M Q ) has 5 C = 20.7 + 
1.9. On the other hand, the relatively low-contrast, third- 
quadrant OCs Ruprecht78 (b — —1.88°; d Q ~ 1.3kpc; 
M - 860 M Q ), Haffner4 (6 = -3.62°; d - 1.7 kpc; M ~ 
260 M Q ) and Trumpler 13 (b = -2.34°; d Q ~ 1.9 kpc; M ~ 
420 M Q ), present 5 C = 4.2 ± 2.6, 6 C = 4.3 ± 2.2, and 5 C = 
8.9 ± 5.5, respectively. 

Completeness is as well affected by cluster distance 
from the Sun in depth-limited photometry. OC stars 
fainter than some apparent magnitude may present ex- 
ceedingly large photometric errors and consequently, end 
up not being separated from the background. To esti- 
mate the dependence of this effect on distance from the 
Sun we consider an artificial OC characterized by a stan- 
dard Salpeter (Salpeter 1955) mass function (0(m) oc 
m~( 1+x \ with x — 1-35) with the turnoff at m = 5Mq. 
Based on our previous experience in the analysis of OCs 
with 2MASS data (e.g. Bonatto & Bica 120051 and refer- 
ences therein) we take stars with J = 15 as representa- 
tive of the average faint-magnitude limit which still can 
be used to derive OC structural parameters. Using the 
mass-luminosity relation derived from the lOMyr, Solar- 
metallicity Padova isochrone (Girardi et al. 120021) . the 
low-mass limit corresponding to J = 15 can be expressed 
as a function of distance from the Sun as mi ow (do) = 
-1.18 + 0.61og(d Q ) + 3.4 x 10- 6 d Q 3/2 . Integration of the 
mass function from the low-mass limit to the turnoff yields 
an estimate of the effective number of stars above the 
background, N*. To a good approximation, this number 
can be expressed as a function of distance from the Sun 
(for 0.2 < d Q (kpc) < 6) asN*(d ) oc (do/lkpc)- 1 . Thus, 
a given OC placed twice as distant from the Sun ends up 
presenting only about half of the stars above the back- 
ground. Roughly speaking, this effect is similar to dou- 
bling the value of at g in eq. ^ which in turn produces a 
smaller value of 6 C . 

From the observational point of view, recent proper- 
motion and parallax data have been used to explore the 
low-contrast end (even less concentrated than Trumpler 
type IV) of the OC distribution. An example is the use of 
HIPPARCOS data to search for poorly-populated OCs or 
candidates in the solar- vicinity (d© < 500 pc) by Platais, 
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Kozhurina-Platais & van Leeuwcn (1998) which yielded 
about 10 such objects which have large angular sizes. One 
must be aware of the completeness effects discussed above, 
but it is also interesting to point out that at least for the 
Solar-vicinity the number of extremely-low contrast OCs 
appears not to be exceedingly large. 

We conclude from the above that the present sample 
(and the results derived therefrom) - although certainly 
affected by completeness - can be taken as statistically 
representative of OCs of Trumpler types I to III. 

The discussions above raise important points related to 
sample completeness and make it clear the need to under- 
take efforts to survey more thoroughly with deeper pho- 
tometry - at least the Solar neighbourhood - to complete 
the local census of the OCs, particularly the low-contrast 
ones. However, they alone do not allow us to objectively 
quantify the completeness either of a particular OC or the 
average completeness of the whole sample. We will return 
to this issue in Sect. El where the quantitative effects of 
completeness will be derived by means of simulation. 

4. The observed disk scale height 

We calculated the vertical position z of all OCs with 
Galactic longitude and latitude, and distance from the 
Sun. The resulting distributions are shown in Fig. for 
the different age groups. The uncertainties correspond 
to 1(T Poisson errors. Similarly to the Galactic latitude 
(Fig. EJ, the z-distributions also depend on OC age, in 
the sense that older OCs are more evenly distributed 
in z than the young OCs. The young (panel (b)) and 
moderate-age (panel (c)) OCs have distributions follow- 
ing exponential-decay profiles, while the old OCs (panel 
(d)) are almost uniformly distributed in height. The total 
number of OCs in each distribution is 402 (young), 148 
(moderate-age) and 86 (old). Note that the number of 
clusters with | z |> 400 pc, which are not shown in Fig.0 
is for the young, 12 for the moderate-age and 27 for the 
old OCs. 

Thus, the near uniformity in the z-distribution of the 
old OCs cannot be accounted for by small-number statis- 
tics alone, but it appears to reflect their intrinsic distribu- 
tion. 

In order to reach a higher z-resolution we built the all- 
ages distribution (panel (a)) with bins of 5pc, while for 
the remaining distributions the z-bin is 10 pc. The higher 
z-resolution is necessary in order to derive more accu- 
rately the distance of the Sun from the Galactic plane 
while at the same time keeping the la Poisson errors non- 
prohibitively large. 

We fit the all-ages distribution in panel (a) with the 
function <f>(z) = </>oe _ " z+z °'/ Zh ' , where z is the displace- 
ment of the Sun with respect to the Galactic plane and 
Zh is the vertical scale height. For this particular analysis 
we center the coordinate system on the Sun. We obtain 
z = 14.8 ± 2.3 pc (above the plane), which within the un- 
certainties, is in excellent agreement with the values cal- 
culated by Cohen (1995) and Hammcrslcy et al. (1995) 




-400 -300 -200 -100 100 200 300 400 



z (pc) 

Fig. 7. Distribution of OCs in terms of height over the 
Galactic plane. Age groups are considered separately. 
Solid lines: least-squares fit with c/)(z) = o e _ ^ z+z °^ Zh ' . 
Shaded areas: la-standard deviation. Top panel: bins of 
5pc in z. Bottom panels: 10 pc bins. 

using stars. This value of z was kept fixed when fit- 
ting the young and moderate-age OCs distributions. The 
results are in Table ^ It is remarkable that large sam- 
ples of different objects (stars and OCs) provide essen- 
tially the same value of the Solar displacement from the 
plane. For completeness we also applied the unconstrained 
exponential-decay fit to the young and moderate-age OCs. 
For the young OCs we obtain z Q = 14.2 ± 2.3 pc and 
Zh = 47.9 ± 2.8 pc which, within the uncertainties are in 
excellent agreement with the constrained values (Tabled . 
However, the unconstrained fit of the less-uniform distri- 
bution of the moderate-age OCs produced a large un- 
certainty in zq, 9.8 ± 11.4 pc and a similar value of Zh, 
145 ± 25 pc. We conclude that the value of the displace- 
ment of the Sun over the Galactic plane derived from the 
distribution of the young OCs is, within the uncertainties, 
the same as the average value which considers OCs of all 
ages. 

According to the above, the young OCs are distributed 
close to the disk with a scale height Zh ~ 48 pc. The scale 
height derived from the moderate-age OCs is ~ 3 times 
larger than that of the young OCs. Our value of Zh for the 
young OCs agrees with that derived by Janes & Phelps 
(1994) for the similar age group. On the other hand, in the 
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Table 1. Observed disk-structural parameters 



Age range 


00 


20 


Zh 


CC 


(Myr) 


/ — i \ 

(PC ) 


(PC) 


(PC) 




(1) 


(2) 


(3) 


(4) 


(5) 


All ages 


5.7 ±0.3 


14.8 ± 2.4 


57.2 ±2.8 


0.87 


< 200 


4.0 ±0.3 


14.8 ± 2.4 t 


47.9 ±2.8 


0.93 


200 - 1000 


0.5 ±0.1 


14.8±2.4 t 


149.8 ± 26.3 


0.64 



Table Notes. 0(z) = O e _l(z+z © )/Zhl . Col.5: correlation coeffi- 



cient, (t): parameter kept fixed. 
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z (pc) 

Fig. 8. z-distribution function of OCs (younger than 
1 Gyr) more internal (panel (a)) and more external (panel 
(b)) than the Solar circle. 



present case the old (age > 1 Gyr) clusters are uniformly 
distributed in z, in contrast to the Zh w 375 pc distribution 
derived by Janes & Phelps l|1994|l . 

4.1. Dependence ofz^ on Galactocentric distance 

In recent years observational evidence in favour of disk 
thickening with increasing Galactocentric distance has 
been presented (e.g. Kent, Dame & Fazio HTjfjTl Janes & 
Phelps 1994). To examine this issue we build z-distribution 
functions for clusters more internal and more external 
than the Solar circle. We exclude from this analysis the old 



OCs because of their nearly-uniform z-distribution (panel 
(d) in Fig. 0|. The results are in Fig. El 

From the fits we conclude that the average scale height 
for disk-regions outside the Solar circle is Zh = 78.1 ± 
5.9 pc (correlation coefficient CC = 0.87), while for regions 
interior to the Solar circle Zh = 39.3 ± 3.3 pc (CC = 0.94). 
This represents an increase in the average Zh by a factor 
of w 2 along the observed Galactocentric distance range. 

5. Sample completeness simulation 

As already suggested by the histograms in Fig. com- 
pleteness effects are most likely affecting critically the 
observed radial distribution of OCs (</>(r)) for large dis- 
tances from the Sun, especially towards the Galactic cen- 
ter/bulge. The discussions in Sect. E3 pointed (i) the crit- 
ical role played by the background stellar density in the 
detection of low-contrast OCs, particularly in directions 
intercepting the Galactic center/bulge, and (ii) the sharp 
decrease of observable stars above the background with 
increasing distance from the Sun. Consequently, the de- 
pendence of completeness on Galactocentric distance (as 
well as on height over the plane) has to be taken into ac- 
count so that more objective results on the intrinsic spatial 
distribution of OCs can be drawn. 

To minimize the effects of distance from the Sun and 
background density we restricted the following analysis to 
the OCs located in the shaded region shown in panel (d) of 
Fig. 2] This region contains the OCs with Galactocentric 
distance in the range 6.7 to 9.3 kpc which on average are 
1.3 kpc distant from the Sun. Considering this distance 
along the Solar circle in both directions the restricted re- 
gion is swept out by azimuthal angles in the range —9.23° 
to ±9.23° with respect to the Galactic center. As a con- 
sequence of the resulting geometry of the region, 90% of 
the OCs included in the restricted sample have distances 
from the Sun of at most d© = 1.3 kpc, and only 3% have 
d in the range 1.5 to 1.8 kpc. 

In Fig. [5] we show the distribution function (0(r) = 
^f^) of OCs (in the restricted sample) in terms of 
Galactocentric distance. As already noted in Fig.^the ob- 
served OC distribution presents a maximum at the Sun's 
position and falls off both for regions internal and exter- 
nal to the Solar circle, with a V-shaped dip at doc ~ 
7.5 — 7.8 kpc (the region before the Sgr-Car Arm), van 
den Bergh & McClurc ( 1980) attributed the sharp drop in 
the number of OCs for regions internal to the Solar circle 
to enhanced OC-disruption through the frequent collisions 
with molecular clouds in the inner regions of the Galaxy, as 
well as dissolution through tidal effects due to the Galaxy. 
However, when compared to the expected distribution of 
an exponential disk, </>(r) oc re~W RD ' (Sect.0, the drop 
in the observed distribution for dec < 8 kpc turns out ex- 
ceedingly large to be accounted for by disruption related 
to dynamical effects alone. 

Considering the above facts we derive the spatial de- 
pendence of completeness by incorporating the qualitative 
arguments presented in Sect. 01 to measurements of the 
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Fig. 9. Observed distribution function of OCs in terms 
of Galactocentric distance in the restricted zone. Dashed 
line: exponential-disk distribution 0(r) oc re~' r / RD '. Note 
the depletion in the number of OCs for doc < 7.8 kpc, 
and the V-shaped dip at dec ~ 7.5 — 7.8 kpc. 

background stellar density as seen by 2MASS photometry 
in different directions around the Sun. 

As a first step we measured the background stellar den- 
sity of the Galaxy as a function of I and b using 2MASS 
photometry. We counted the number of stars brighter than 
J = f 5 contained in circular areas with f ° in radius for the 
Galactic latitudes b = 0°, ±2.5°, ±5°, ±10°, ±25° and ± 
50° in steps of 5° in I. The selected circular area is large 
enough to accommodate most of the local fluctuations in 
the number of stars. We illustrate some of the northern 
background density profiles in panel (a) of Fig. where 
the marked difference in the average background in direc- 
tions towards the Galactic center/bulge with respect to 
those towards the anti-center is clear. The corresponding 
southern profiles are very similar to the northern ones. 
The background contribution for a given OC (i.e. at a 
specific pair of the coordinates I and b) is obtained by 
interpolation among the above curves. 

Our analysis is based on the detectability of OCs - par- 
ticularly the more distant ones - projected towards differ- 
ent regions of the Galaxy, a process which depends essen- 
tially on the density contrast of the central region of a 
cluster with respect to the field, fn this sense we assume 
a typical angular radius of 1' to characterize the central 



region of the clusters. To estimate the effective number 
of stars of the OCs in the restricted zone we counted the 
number of stars brighter than J = 15 (~ er(0)) within 
the R = 1' central region, regardless of the OC distance 
from the Sun. The corresponding background contribu- 
tion is subtracted therefrom to yield the effective number 
of stars, N* = <r(0) — Obg = &ok- 

To avoid the high background density towards the 
Galactic center/bulge we restricted this procedure to the 
OCs external to the Solar circle. We assumed the WEBDA 
coordinates as the center of the OCs from which we ex- 
tracted the photometry. The resulting histogram of the 
number of OCs with a given effective number of stars is 
shown in panel (b) of Fig. (solid line) . For a more ob- 
jective comparison of the effective stellar content among 
the OCs we recalculated N* assuming d© = 500 pc for 
all OCs, using the relation N * (d Q ) ~ d© _1 (Sect. |3J). 
The resulting histogram is in panel (b) of Fig. EH (dotted 
line). Interestingly, the corresponding distribution func- 
tion (panel (c)) is well represented by an exponential- 
decay profile, </> 50 o(N*) oc e -( N */ 22 ). 

OCs with off-center coordinates in WEBDA and those 
with large projected areas such as the Hyades and Pleiades 
would simply enhance the low-N* tail of the </>5oo(N*) dis- 
tribution (panel (c) of Fig. I10|l. which would not change 
appreciably the shape of the distribution. 

Our basic assumption is that the observed distribu- 
tion of the number of OCs with a given effective number 
of stars outside the Solar circle (</>soo(N*)) both represents 
the intrinsic distribution of OCs in the Galaxy and is es- 
sentially isotropic with respect to the Sun. Based on this 
we randomly assign to each observed OC position (£, b, 
d©) in the restricted zone a value of N* (corresponding to 
d Q = 500 pc) with a number-frequency distributed accord- 
ing to the observed histogram (panel (b) of Fig. fTTTft - and 
recalculate N* for the actual distance from the Sun of the 
OC. Finally, we divide this last number by the correspond- 
ing number of background stars (calculated at the respec- 
tive I, b) to obtain the contrast parameter 5 C = 1 + ^^ 3 - 
We count as detectable only the OCs with 6 C > 2 (i.e. 

> 1). Completeness at a given position is then the 
fraction of retrieved OCs with respect to the input ones. 
To test convergence of this procedure we ran 10 5 simula- 
tions. We illustrate the convergence of this process in panel 
(d) of Fig. ITUI where we plot histograms of the resulting in- 
tegrated completeness value with respect to the number of 
simulations, both for the OCs inside and outside the Solar 
circle. The histograms are well represented by gaussian 
distributions centered at 31.7% and 63.8%, with standard 
deviations of 2.8% and 3.5%, respectively for the OCs in- 
side and outside the Solar circle. As expected because of 
the higher background density, the average completeness 
for the OCs inside the Solar circle is about half of that for 
the OCs outside it. The average integrated completeness 
throughout the entire restricted zone is 47.3 ± 2.1%. 
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Fig. 10. Panel (a): average background density (2MASS 
stars brighter than J = 15) to the north of the Galactic 
plane. la-Poisson errors are smaller than the symbols. 
Panel (b): histograms of the number of OCs (outside the 
Solar circle in the restricted zone) with a given effective 
number of stars at the actual distance from the Sun (solid 
line) and at d© = 500 pc (dotted line). Panel (c): the 
distribution function at d© = 500 pc can be fitted by 
05Oo(N*) oc e~' N */ 22 ); Shaded area: la standard devia- 
tion from the fit. Panel (d): average integrated complete- 
ness for the OCs inside (solid line) and outside (dotted 
line) the Solar circle and throughout the restricted zone 
(dashed line). 



The completeness fraction as a function of 
Galactocentric distance (for the OCs in the re- 
stricted zone) is shown in panel (a) of Fig. ^] As 
expected, the completeness is highest around the Sun 
(fcomp(8kpc) « 0.68) and drops off towards external 
regions and more sharply towards the Galactic center. 
Dividing the observed OC distribution function (Fig. |UJ) 
by fcomp(r) produces a completeness-corrected radial 
distribution which follows the expected exponential- 
disk profile (panel (b) of Fig. 111(1 along the entire 
Galactocentric distance range of the restricted zone, 
with a scale-length Rd = 1.4 ± 0.2 kpc. This estimate 
of Rd represents about 40% of the stellar overall disk 
scale-length (de Vaucouleurs & Pence I1978fl . Our de- 
termination probably reflects the Solar-vicinity value of 
Rd and/or that the OCs distribute in the disk with a 
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Fig. 11. Panel (a): completeness fraction as a function 
of Galactocentric distance in the restricted zone. Panel 
(b): observed (empty circles) and completeness-corrected 
(filled circles) OC distributions. Dashed-line: fit with an 
exponential-disk with a scale-length Rd = 1.4 ± 0.2 kpc. 



different scale length than the stars. Remark that the 
V-shaped dip at dec ~ 7.5 — 7.8 kpc (attributed to the 
edge before the Sgr-Car Arm - Sect. I2.1f> is still present 
in the completeness-corrected distribution function. In 
this sense, the above completeness-simulation procedure 
is sensitive as well to conspicuous physical structures of 
the Galaxy. 

Since our analysis is based on magnitude-limited 
2MASS photometry (J = 15), excess reddening of back- 
ground sources towards the inner Galaxy with respect 
to those in the anti-center (Fig. [SJ may have underes- 
timated the contrast parameter of the OCs in that re- 
gion because of the increased differential brightness of the 
brightest stars. As a consequence, completeness for the 
OCs in the inner Galaxy would be underestimated, which 
might account for the over-correction of 0(dcc) with re- 
spect to the exponential-disk fit, as suggested by Fig. ^2 
for doc < 7.5 kpc. 

Considering the completeness-corrected OC distribu- 
tion, the fraction of observed OCs in the restricted zone 
corresponds to about 47% of the probable total (essen- 
tially of Trumpler types I to III - Sect. 
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6. Extension to the whole OC sample 

In the previous section we corrected for completeness 
the observed radial distribution of OCs in the restricted 
zone and, as a result obtained the expected exponential- 
disk profile. At this point it is interesting to check 
whether the completeness-simulation procedure retrieves 
an exponential-disk profile when applied to the whole OC 
sample. We test this hypothesis considering 2 different 
simulation approaches. 

6.1. Actual-position simulation 

In this simulation the effective number of stars (Sect. |5j 
N* is randomly distributed (with the number-frequency 
according to the histogram in panel (b) of Fig. I10|l among 
the 654 OCs assuming the WEBDA spatial locations (£, 
b, d©). The completeness fraction of each OC is then the 
average value of the retrieval rate after running 10 5 sim- 
ulations. Similarly to the restricted-zone simulation, the 
individual ler-standard deviation of the completeness frac- 
tion is at most 3% of the corresponding average value. 

The observed radial distribution of OCs in the range 
5 < dec (kpc) < 14 (Fig. 112(1 presents a maximum 
at the Solar circle and drops off both towards larger 
Galactocentric distances and more sharply to the Galactic 
center. The radial completeness fraction resulting from 
the actual-position simulation (Fig. 113(1 presents simi- 
lar features. Correction of the observed distribution for 
completeness produces a radial distribution which follows 
an exponential-disk profile (panel (a)) characterized by a 
scale-length of Rd = 1.5 ± 0.1 kpc (correlation coefficient 
CC = 0.88). 

The average integrated completeness throughout the 
disk is 33.2 ± 1.5%. 

6.2. Random-position simulation 

In this case we build 1000 artificial test positions (£, b, 
d©) over which we randomly distribute the effective num- 
ber of stars as in the actual-position simulation. The £ 
coordinate is randomly selected around the Sun while b 
is drawn according to the observed number-frequency his- 
togram shown in panel (e) of Fig. |31 to match most of the 
observed distances from the Sun (panel (a) of Fig. 0J, d Q 
is randomly selected from the range - 7 kpc. Because of 
the random nature of the test OC positions, the resulting 
radial completeness fraction (Fig. 113(1 turns out similar 
to, but considerably smoother than that derived in the 
actual-position simulation. 

The completeness-corrected radial distribution of OCs 
(panel (b) of Fig. 112(1 follows an exponential-disk pro- 
file characterized by a scale- length Rd = 1.9 ± 0.1 kpc 
(cc = 0.88), about 30% larger than that obtained with 
the actual-position simulation, within the uncertainties. 
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Fig. 12. Completeness-correction of the observed radial 
distribution of OCs (empty circles) according to the 
actual-position (top panel) and random-position (bot- 
tom panel) simulations. The resulting radial distributions 
(filled circles) follow exponential-disk profiles character- 
ized by scale-lengths Rd ~ 1.5 kpc and Rd ~ 1.9 kpc, 
respectively. 

6.3. Comparison of simulations 

The simulations discussed above are based on concep- 
tually different approaches with respect to the spatial 
distribution of OCs. Despite this, they produced sim- 
ilar completeness-corrected radial distributions of OCs, 
which enhances the significant role played by complete- 
ness in depth-limited photometric surveys throughout the 
Galaxy. However, the completeness-corrected profile ob- 
tained with the random-position simulation produces a vi- 
sually better overall fit with the expected exponential-disk 
profile fFig. 112(1 as compared to the actual-position simula- 
tion, especially near the inner and outer dec borders. The 
completeness-corrected distributions derived from both 
simulation approaches present similar disk scale lengths, 
Rd ~ 1.5 — 1.9 kpc, about 43% - 54% of the canonical 
value (Rd ~ 3.5 kpc) derived with stars by de Vaucouleurs 
& Pence (1978). This fact may reflect intrinsic differences 
in the spatial distribution of OCs with respect to stars. 

In Fig. we compare the radial completeness frac- 
tions resulting from both simulations. They basically 
agree, except for dec < 6.6 kpc where the actual-position 
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Fig. 13. The actual-position radial completeness fraction 
(solid line) drops off more sharply towards the Galactic 
center for dec < 6.6 kpc than that of the random-position 
simulation (dot-dashed line). Inset: the distribution of 
the number of OCs with a given effective (background- 
subtracted) number of stars for the whole OC sample is 
well fitted by an exponential-decay profile. Shaded area: 
la standard deviation from the fit. 
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Fig. 14. Completeness fraction as a function of z. Left 
panel: completeness fraction for OCs internal (dashed line) 
and external (dotted line) to the Solar circle, and the aver- 
age disk completeness (solid line). Right panel: the shape 
of the average disk completeness results from the com- 
bination of OCs at different ranges of distance from the 
Sun. 



completeness fraction drops off more sharply towards the 
Galactic center than that from the random-position simu- 
lation. Considering that in the random-position simulation 
the OCs are uniformly distributed throughout the Galaxy, 
as opposed to observed OC coordinates in the actual- 
position, perhaps this difference may suggest that the 
tidal-disruption effects due to the Galactic center/bulge 
begin to be observationally significant at dec ~ 6.6 kpc. 

Finally, in the inset of Fig. ^| we show the distribution 
of the number of OCs with a given effective (background- 
subtracted) number of stars for the whole OC sample. 
Similarly to the OCs in the restricted zone this distribu- 
tion follows an exponential-decay profile, but with a char- 
acteristic number of stars of AN* « 95. The decay-scale 
for the whole sample is about 4 times larger than that de- 
rived for the OCs outside the Solar circle in the restricted 
zone (Sect. |SJ| . This discrepancy can be accounted for by 
an observational bias, since detection probability in direc- 
tions intercepting the Galactic center/bulge is higher for 
the more populous OCs. Interestingly, the analytical ex- 



pression of the overall ^500 resembles that of the galaxy 
luminosity function f Schecht er 1 1 9 76[1 . 

6.4. Completeness in vertical distributions 

In Sect. I4.ll we showed that the observed scale height de- 
pends on Galactocentric distance, in the sense that the OC 
distribution outside the Solar circle is about 2x broader 
than that inside it. At this point it may be interesting to 
determine how completeness affects the vertical distribu- 
tion and derive its intrinsic shape. 

In Fig. E| (left panel) we show how the completeness 
fraction as a function of z depends on Galactocentric dis- 
tance, centered at the Solar position over the disk (z Q = 
14.81 pc, Sect.0}. The curves have been produced with the 
random-position simulation fSect. l6~2)l . As expected, the 
OCs internal to the Solar circle suffer more severely from 
completeness than the external ones. Interestingly, both 
curves present similar shapes, with a peak at z = z Q , sharp 
declines for increasing |z| with a minimum at |z| « 340 pc, 
and an increase for larger values of |z|. As shown in the 



C. Bonatto et al.: Spatial distribution of open clusters 



13 




z(pc) 

Fig. 15. Completeness-corrected z-distributions. Panel 
(a): overall radial distribution (scale height Zh ~ 78 pc). 
Panel (b): OCs internal to the Solar circle (zj, w 85 pc). 
Panel (c): OCs external to the Solar circle (zh ~ 119 pc). 
Empty circles: observed OCs. Filled circles: completeness- 
corrected z-distribution. 

right panel OCs at different ranges of distance from the 
Sun present similar dependence with z, with a minimum at 
z = zq and an increase for larger |z| . However, the absolute 
values of completeness drop for more distant OCs. The fi- 
nal shape of the completeness fraction results from the 
combination of OCs at different distances from the Sun. 
In particular, the peak at z = z Q can be accounted for 
by the high completeness of the OCs closer than w 1 kpc 
from the Sun. 

The observed z-distributions of OCs internal and 
external to the Solar circle, as well as the over- 
all one have been corrected for completeness with 
the corresponding curves derived above. The resulting 
completeness-corrected distributions along with the re- 
spective exponential-decay fits are shown in Fig. 1151 The 
completeness-corrected scale heights are Zh = 78 ± 3 pc 
(observed: Zh w 57 pc) with correlation coefficient CC = 
0.88, z h = 85 ± 13 pc (observed: z h « 39 pc) with CC = 
0.77, and Zh = 119 ± 12 pc (observed: Zh ~ 78 pc) with 
CC = 0.75, respectively for the overall, internal and ex- 
ternal to the Solar circle distributions. Consequently, com- 
pleteness correction of the z-distributions preserved the 
dependence of scale height on Galactocentric distance, in- 



creased the absolute values of zj, with respect to the ob- 
served ones and decreased the observed ratio of Zh out- 
side/inside the Solar circle from 2 to 1.4. 

Completeness-correction of the observed z- 
distributions of the OCs younger than 200 Myr and 
those with age in the range 200 Myr - 1 Gyr (panels 
(b) and (d) of Fig. |7| respectively) increased the scale 
height by a factor of ~ 1.42. Consequently, the young 
OCs distribute vertically following an exponential-decay 
profile with scale height Zh ~ 68 pc. For the 200 Myr 

- 1 Gyr OCs the completeness-corrected scale height 
is Zh ~ 230 pc. Considering the observational uncer- 
tainties, these values are comparable to the Hi scale 
height in the Solar neighbourhood, which we take as 

- ±FWHM(Hi) w 115pc (Dickey & Lockman ITMTjl . 
This must reflect the association of star formation and 
the parent gas. The same conclusion applies to CO, for 
which surveys of the outer Galaxy (Heyer et al. ^998), 
where the Perseus arm lies, find an overall mean CO scale 
height of 113 pc. 

6.5. Age distribution function 

The histogram with the number of observed OCs in bins 
of age (Fig. EJ was converted into the observed age- 
distribution function (panel (a) of Fig. 1161) . This function 
can be fitted by a combination of exponential-decay pro- 
files with age scales of T young = 123 ± 13 Myr (characteris- 
tic of the young OC population) and r id = 2.4±1 Gyr (old 
OCs). The correlation coefficient of the fit is CC = 0.94. 
Clearly the observed old OC population represents an 
excess over the extrapolation to old ages of the young 
OCs distribution function. An additional excess in the ob- 
served age-distribution function shows up in the age range 
~ 600 Myr to ~ 1.5 Gyr, which can probably be identified 
with the moderate-age OC population. The present young 
and old age time-scales are comparable to those derived by 
Janes & Phelps (.1994), 200 Myr and 4 Gyr for the young 
and old OCs, respectively. 

The observed age-distribution function (panel (a) of 
Fig. I16fl basically contains the survivor clusters after 
emerging from the parent molecular clouds. Lada & Lada 
(2003) estimate that, for embedded clusters within 2 kpc 
from the Sun, only a fraction of 4-7% emerge as OCs. Bica, 
Dutra & Barbuy (2003) compiled 276 embedded clusters 
(mostly infrared clusters) and stellar groups from the lit- 
erature located mostly within 3 kpc from the Sun. From 
the age-distribution function in Fig. 1161 we estimate that 
the number of observed OCs with ages in the ranges 5- 
10 Myr, 10-15 Myr, and 15-20 Myr are 21 ±3, 21 ±3 and 
20 ± 3, respectively. Considering an age range of 5 Myr 
for the embedded clusters (typical lifetime of H II regions) 
and comparing with the number of OCs in the subsequent 
5 Myr-bins (see above) we can derive a survival rate of 
~ 8%, similar to Lada & Lada I20U3| 's upper limit. More 
recently Dutra et al. ( 2003) and Bica et al. ( 2003 ) surveyed 
for new embedded clusters, groups and candidates in the 
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Fig. 16. Panel (a): the observed age-distribution can be 
fitted with two exponential-decay laws with time-scales of 
r youn g = 123 ± 13 Myr (dotted line) and T id = 2.4 ± 1 Gyr 
(dashed line) characterizing the young and old OCs, re- 
spectively. Panel (b): same as (a) for the completeness- 
corrected age distribution; time-scales are Ty OU ng = 89 ± 
11 Myr (dotted line) and r id = 1.3 ± 0.5 Gyr (dashed 
line). Shaded region: ler-standard deviation of the com- 
bined young+old fit. 

direction of radio and optical H 11 regions and molecular 
clouds. These objects are typically within ~ 5 kpc from the 
Sun according to kinematical distances of the related H 11 
regions. Adding these newly discovered embedded clusters 
to those of Bica, Dutra & Barbuy l|200 3). totaling 622 
objects, the survival rate drops to ~ 3.4%. The present 
estimates match those of Lada & Lada {2003 ) . 

One output of the actual-position simulation 
(Sect. EJ is the average individual completeness of 
each OC in the present sample. The reciprocal of the 
completeness is essentially the detection probability of 
that particular OC - and consequently of its intrinsic 
parameters, such as the age. In this sense we can apply the 
same procedure to correct the observed age distribution 
for completeness. The resulting completeness-corrected 
age distribution (panel (b) of Fig. 1161) retains the basic 
features of the observed one, including the intcrmcdiate- 
age OC excess. However, because the distribution in 
Galactic latitude of the young OCs is narrower than that 
of the old OCs (Fig. OJ), the young OCs are subject on 



average to higher incompleteness. Consequently, com- 
pleteness correction should enhance the number of young 
OCs relative to the old ones. Indeed, the observed ratio 
of the number of OCs younger than 1 Gyr with respect to 
the older ones is « 7, while in the completeness-corrected 
distribution this value increases to ~ 14. 

Similarly to the observed profile, the best-fit of the 
completeness-corrected age distribution results from a 
combination of exponential-decay curves with time scales 
Tyoung = 89±11 Myr and t oM = 1.3±0.5 Gyr (CC = 0.94), 
respectively for the young and old OCs. The completeness- 
corrected time scales result comparable to the observed 
ones, within the uncertainties. We note however that the 
resulting time scales of both OC populations correspond to 
about half of the corresponding values estimated by Janes 
& Phelps <|1994f) . This might suggest faster destruction- 
rates than those implied by Janes & Phelps ljl994|) . and 
particularly with respect to the ~ 600 Myr destruction- 
time scale suggested by Bergond, Leon & Guibert H2001fl . 

From the above we conclude that evolutionary scenar- 
ios based on constant rates of cluster formation and dis- 
ruption do not apply to the double exponential-decay age 
distribution. Probably a time- varying star-formation rate, 
such as a massive initial burst at the Gyr scale, might fit 
the data. In this context, the suggestion of a galaxy merger 
in the early phases of the Milky Way (Chen, Stoughton, 
Smith et al. I2001[) and the subsequent disturbances to 
the disk and burst of star formation, could account for 
the observed excess in the number of old OCs surviving 
to the present time. One clue to disentangle the cluster 
formation/disruption-rate problem may be provided by 
the systematic study of OC remnants (e.g. Pavani et al. 



6.6. Galactocentric distance vs. age 

In Fig. El we examine the relation of Galactocentric dis- 
tance with age. We divided the OC sample in bins of 
age and calculated the average age and dec in each 
bin, as well as the respective standard deviations. Lynga 
lP^8"2) and Tadross et al. ijM^I concluded that older 
clusters are found preferentially in the outer parts of 
the Galaxy, while younger clusters are more evenly dis- 
tributed. The observed OC distribution is in panel (a) and 
the completeness-corrected one is in panel (b). The large 
l<7-standard deviation bars in the observed distribution, 
particularly in dec, preclude any definitive conclusion on 
this issue. 

However there appears to exist a trend in which older 
clusters are found preferentially at larger Galactocentric 
distances. This trend is enhanced after applying the com- 
pleteness correction to the observed distribution (panel 
(b)). A linear-least squares fit to the completeness- 
corrected distribution results in the relation (dec) = 
(7.83 ± 0.27) + (0.0012 ± 0.0003)(Age) (CC = 0.87), for 
doc m kpc and Age in Myr. 
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Fig. 17. Relation of Galactocentric distance and age. 
Panel (a): observed OCs. Panel (b): completeness- 
corrected. The bars show the lc-standard deviations 
of the average values of age and Galactocentric dis- 
tance within the respective bins. Dashed line: linear-least 
squares fit to the completeness-corrected points. Shaded 
region: lcr-standard deviation of the fit. 

7. Estimate of the total number of open clusters 

We describe the spatial number-density of open clusters 
in the disk as a combination of exponential-decay laws for 
the radial and vertical components, p(r, z) — dN ° c — 



Po e 



-(r/R D ) p — |z/z h [ 



rdrdz d6» 



, where Rd is the disk scale-length, Zh 
is the scale-height and p is the number-density of OCs 
at the Galactic center. Using cylindrical coordinates and 
defining D as the Galactocentric distance of the Sun, 
z the displacement of the Sun from the Galactic plane, 
and pq the local (solar vicinity) number-density of OCs, 

we can write p(r, z) = p Q e v r d / e V h / . Integration 
over the z-axis and through a region with opening angle 
Af? produces the radial distribution function of OCs 



dN ol 
dr 



4-7T 
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p Q r z h e z h e 
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The total number of open clusters in the disk is then 
(A0 = 2tt) 



fa (Ee) o 
47rpoZh e z h e\ RD 'R D 



1 



1 



Rd 



(3) 



where Rg is the disk radius (Rg ~ 15kpc, Binney & 
Tremaine 1987) 4 . Accordingly, estimates of both the solar- 
vicinity OC density and the local disk scale length can be 
obtained by fitting eq. 0to completeness-corrected distri- 
bution functions. 

Based on eq. F2] and the completeness-corrected dis- 
tribution functions of the actual-position and random- 
position simulations for the whole OC sample (with A£ ~ 
0.11, Fig. we derive p Q = 759 ± 34 kpc -3 and p Q = 
774 ± 33kpc~ 3 , respectively (Fig. El) . 

Extrapolation of the above values of po and Rd to 
the whole disk using eq. produces a total number of 
OCs in the Galaxy of - 3.7 x 10 5 and - 1.8 x 10 5 , 
respectively for the actual-position and random-position 
completeness-corrected distribution functions. These es- 
timates were obtained by direct extrapolation of the 
completeness-corrected distributions down to the Galactic 
center. Consequently, they do not take into account de- 
pletion in the number of OCs by tidal disruption and en- 
hanced frequency of collisions with molecular clouds in the 
inner parts of the Galaxy. In this sense, the above numbers 
of Galactic OCs must be taken as upper-limits. 

For the completeness-corrected distribution function in 
the restricted zone (with a scale-height of 57 pc (Table 
and |f = 0.05) we derive R D = 1.4 ± 0.2 kpc and p Q = 
795 ± 70 kpc -3 fFig. Illfl . This number can be compared 
with that obtained assuming a uniform, cylindrically- 
symmetric distribution of OCs around the Sun. With a 
total number of OCs up to a distance from the Sun of 
1.3 kpc of 341, and using a vertical height of as 2x the 
scale-height we derive p© = 281 ± 16 kpc -3 , which cor- 
responds to the observed OCs. Correcting for the aver- 
age completeness in the restricted zone (~ 47%) produces 
a probable total of ~ 730 OCs, and a number density 
of pq = 600 ± 33 kpc -3 , which within the uncertainties 
roughly agrees with the value derived from the distribu- 
tion function. 



8. Concluding remarks 

In the present paper we use a sample of 654 open clus- 
ters (typically Trumpler types I to III) with published 
parameters to derive statistical properties related to age, 
distance from the Sun, reddening, Galactocentric distance 
and height with respect to the Galactic plane. 

The population of OCs younger than 200 Myr dis- 
tributes vertically in the disk following an exponential- 
decay profile with a scale height of Zh = 47.9 ± 2.8 pc. 
Clusters with ages in the range 200 Myr to 1 Gyr dis- 
tribute vertically with Zh = 150 ± 27 pc, while older clus- 
ter distribute nearly uniformly in height from the plane 
so that no scale height can be derived from exponential 
fits. The average scale height, considering clusters of all 
ages, is Zh = 57.2 ± 2.8 pc. The above OC scale heights 



Because of the exponential factor in eq. 3 e V R o J with 
Rd = 1.5 — 1.9 kpc fSects. UTTl and HOI . our estimates of N oc 
are essentially insensitive to values of Rg > 10 kpc. 
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are considerably smaller than those attributed to the thin 
disk (zh « 0.6 kpc) as derived by means of stars (e.g. de 
Boer et al. 119971 Altmann et al. 120041 Kaempf, de Boer 
& Altmann 2005]). We confirm previous findings that z^ 
increases with Galactocentric distance, being about twice 
as large in regions outside the Solar circle than inside it. 

The asymmetry in the vertical distribution of OCs al- 
lowed us to derive the displacement of the Sun above the 
Galactic plane as zq = 14.8 ± 2.4 pc, which agrees with 
previous determinations using stars. 

Sample completeness affects critically the detection 
particularly of low-contrast OCs internal to the Solar cir- 
cle and those distant from the Sun. We simulate the ef- 
fects of completeness in all directions and for any dis- 
tance from the Sun assuming that the intrinsic distri- 
bution of the number of OCs with a given number of 
stars (above the background) measured in a restricted 
zone outside the Solar circle is isotropic. As a conse- 
quence we derived the radial dependence of complete- 
ness. The observed radial distribution of OCs in terms of 
Galactocentric distance does not follow the expected expo- 
nential profile, instead it falls off both for regions external 
to the Solar circle and more sharply towards the Galactic 
center. Correction for completeness produced radial distri- 
butions which agree with exponential disks throughout the 
Galactocentric distance range 5-14 kpc, with scale lengths 
Rd = 1.5 — 1.9 kpc. These values of Rd are smaller than 
those implied by stars, which may reflect intrinsic differ- 
ences in the spatial distribution of OCs and stars. 

With respect to the vertical disk structure described 
by the OCs, we showed that older clusters distribute in ex- 
ponential disks with larger scale heights than young OCs. 
In addition we confirmed previous observations that the 
scale height increases with Galactocentric distance. The 
latter result was further enhanced by the completeness- 
corrected z-distributions. The completeness-corrected OC 
scale height is comparable to those implied both by the 
Hi and CO distributions. 

We derived a number-density of solar-neighbourhood 
(with distances from the Sun d© < 1.3 kpc) OCs of Pq = 
795 ±70 kpc -3 , which implies a total number of (basically 
Trumpler types I to III) OCs of - 730 of which - 47% 
would already have been observed. Inferences on the to- 
tal number of OCs in the disk can be made based on the 
local number density and radial/azimuthal completeness- 
corrected OC distributions. This estimate is important 
for the realization of what to expect in future surveys 
of the Galaxy. Direct extrapolation of the completeness- 
corrected distributions down to the Galactic center sug- 
gests that the total number of OCs in the Galaxy is in the 
range (1.8 — 3.7) x 10 5 . These numbers must be taken as 
upper-limits because the extrapolation does not take into 
account depletion in the number of OCs by tidal disrup- 
tion and enhanced frequency of collisions with molecular 
clouds in the inner parts of the Galaxy. 

The observed and completeness-corrected OC age- 
distribution functions can be described by a combination 
of two exponential-decay profiles characterizing the young 



and old OC populations with age scales of ~ 100 Myr 
and ~ 1.9 Gyr, respectively. As a consequence, scenarios 
based on constant star-formation and OC-disruption rates 
do not apply to the data. Comparing the number of em- 
bedded clusters from IR studies in the literature and the 
present observed age-distribution function of very young 
OCs we derive survival rates in the range 3.4-8%, in agree- 
ment with the estimates of Lada & Lada ( 2003 ) . 
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